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ABSTRACT 

Context. During the past decade circumbinary disks have been discovered around various young binary stars. Hydrodynamical calcu¬ 
lations indicate that the gravitational interaction between the central binary star and the surrounding disk results in global perturbations 
of the disk density profile. 

Aims. We study the observability of characteristic large-scale disk structures resulting from the binary-disk interaction in the case of 
close binary systems. 

Methods. We derived the structure of circumbinary disks from smoothed-particle hydrodynamic simulations. Subsequently, we per¬ 
formed radiative transfer simulations to obtain scattered light and thermal reemission maps. We investigated the infiuence of the binary 
mass ratio, the inclination of the binary orbit relative to the disk midplane, and the eccentricity of the binary orbit on observational 
quantities. 

Results. We find that ALMA will allow tracing asymmetries of the inner edge of the disk and potentially resolving spiral arms if the 
disk is seen face-on. For an edge-on orientation, ALMA will allow detecting perturbations in the disk density distribution through 
asymmetries in the radial brightness profile. Through the asymmetric structure of the disks, areas are formed with a temperature 2.6 
times higher than at the same location in equivalent unperturbed disks. The time-dependent appearance of the density waves and spiral 
arms in the disk affects the total re-emission flux of the object by a few percent. 
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1. Introduction 

During the past three decades the formation of single stars was 
explored in detail (e.g., see the review |McKee & Qstriker|2007] ) 
Circumstellar disks, which are formed by the conservation of an¬ 
gular momentum during the molecular cloud collaps, play a key- 
role in this process (e.g., [Andrews & Williams ||2005[| Williams | 

& Andrews|2006t [Andrews & Williams |2007[|Isella et al.|2009t 

Andrews et al.||200^ 201Q| ) The lifetime of these disks, which 

have been observed from the optical to the millimeter wave 
length range, is limited to a few million years (|Natta et al 


2000l[Mathieu et al.|2000l|Haisch et al.|2QQll|Wolf et al.|2012i r 

In the optical and near-infrared regime, the scattering of stel¬ 
lar radiation on the surface dominates the appearance of these 
disks, while the thermal reemission is stronger at far-infrared to 
millimeter wavelengths. High-angular observations across this 
wavelength range are important for the analysis of the disk struc¬ 
ture and composition as well as for their evolution (e.g.,|Sauter| 
& Wolf|2011i|Grafe et al.|2013D 
^ However, since the mid-1990s it became evident that a large 
portion of stars forms as binaries (e.g ., [Kraus & Hillenbrand| 
2009[ [Bonnell & Bate|p^94{ |Duch6ne |1999| ) Meanwhile, gas 


and dust disks have been discovered in several young binary sys 
terns (e.g., IGuilloteau et al.|199 9t McCabe et al.|2QQ2||Duchene 


et al. 2007 Romero et al. 2012| ) They come in two flavors: First, 
the systems where at least one component of the binary hosts a 
circumstellar disks, as studied by [Regaly et al.| ( |2011a|b| ) ; and 


second, bin aries with a circumbinary disk (e.g ., GG tau |Pietu 
et aipOlT] V4046 Sgr., [Rodriguez et ar]|2Q10| ). In the second 
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case (circumbinary disks), the orbital motion of the central bi¬ 
nary is predicted to induce periodic perturbations, spiral density 


Kocsis et al.l 


Larwood & Papaloizou[[1997[ 

Gunther & Kley[[2002 Gunther 

et al.[2004 Sotnikova & Grinin 

2007 

[Hanawa et aL|2010[ 

Kaig- 

orodov et al.| 

2010 [Alexander[ 

2012 

[de Val-Borro et al.[ 

20TT] 


2012) Characteristic features resulting from this 


interaction are inner cavities with the radius of about twice the 
semi-major axis of the binary (e.g., [Arty mowicz & Lubow|1994t 
Gunther & Kley|2002| ) This binary-induced inner cavity results 


in a flux-deficit in the mid-infrared wavelength range and a large 
inner hole in the brightness distribution of these objects (at least 
at sub-to-millimeter wavelengths). This appearance is similar to 
what has been prediced and found in the case of transitional 
disks (e.g., Artymowicz & Lubow|1994[[Evans et al.|2009t|Cieza 


[et al.[[2010 ) The existence of circu mbinary planets wa s shown 

through the discovery of Kepler 34b ( [Welsh et al.|2012[ ) , for in¬ 
stance. However, the effect of the binary-induced perturbations 
on the process of planet formation in circumbinary disks is still 
under debate (e.g.,[Paardekooper et al.|2Q12t[Lines et al.|2014t 
[Thebault & Haghighipour|2014[ ) Observational constraints that 
allow verifying the physical processes of the predicted binary- 
disk interactions are required. For this purpose, we investigate 
the feasibility of observing characteristic binary-induced pertur¬ 
bations in circumbinary disks. In particular, we focus on the case 
of embedded close binary systems. Previous studies have shown 
that the orbital motion of binary-induced large-scale density per¬ 
turbations in the disks ca n cause a large-amplitude brightness 
modulation of the system ([Sotnikova & Grinin|20()7] [Demidova 
et al.|2010b} [Demidova et al.|201Qa[ [Grinin et al.|2QlQ[ ). In par- 
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Fig. 1. Illustration of the projected disk density distribution in a face-on 
orientation for models with an eccentric binary orbit. The binary mass 
ratio q is shown in the upper left corner of every subplot. The filled 
circles indicate the positions of the binary stars. 


ticular, the perturbation of the inner part of the disk affects the 
illumination of its outer regions, resulting in an asymmetric im¬ 
age of the disk if observed face-on ( [Demidova et al.|2QT3] ) 

Based on selected circumbinary disk models resulting from 
3D smoothed-particle hydrodynamic (SPH) simulations, we per¬ 
form radiative transfer calculations to derive the thermal struc¬ 
ture and observable quantities of these circumbinary disks. To 
trace the densest regions in the disks best - and thus the most 
prominent density perturbations - one has to make sure to ob¬ 
serve in a wavelength range where they are optically thin. Com¬ 
bined with the second requirement, namely to spatially resolve 
structures on the sub-lOOAU scale, we expect the Atacama Large 
Millimetre/submillimeter Array (ALMA) to be the potentially 
best observatory for future studies. For this reason we prepared 
and analyzed simulated observations of circumbinary disks with 
this interferometer. 

The paper is divided into three parts. First, we describe the 
SPH simulations in Sect.|2]and the radiative transfer calculations 
in Sect.j^ All assumptions and the parameter setup are addressed 
in these sections. Predicted ALMA observations are presented in 
Sect.H 

2. Circumbinary disk structure 

As the basis of our study, the periodically changing gravitational 
interaction between the embedded binary system and the sur¬ 
rounding disk was calculated with 3D SPH simulations. In the 
following, we explain the underlying assumptions and parameter 
setup of our study. We applied the SPH code developed by jSot-j 
|nikova|(|1996|)i for a more detaile d de scription of the SPH mod - 
els see [Demidova et aL] ( |201Qa| ) and [Demidova et aL] ( [2010b[ ) 
For the follow-up radiative transfer calculations (see Sect. [^ we 
assumed a dust-to-gas mass ratio of 1 : 100. 

We focus on a simple disk setup and selected orbital elements 
of the binary and the binary mass ratio. 

Central binary: We assumed a close binary with a semi-major 
axis of a = 2.0 AU, which corresponds to a period of 1.7 yr 
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Fig. 2. Illustration of the projected disk density distribution in a face-on 
orientation for models with a binary orbit that is inclined relative to the 
initial disk midplane. The binary mass ratio q is shown in the upper left 
comer of every subplot. The filled circles indicate the positions of the 
binary stars. 


(^prim — 2.5 Mq). 

Circumbinary disk: We assumed a fully isothermal disk. The ef¬ 
fective viscosity of the disk ( [Shakura & Sunyaev|1973[ ) , was set 
to = 0.03 (for details of implementating the viscosity in our 
simulations see [Sotnikova & Grinin[ ( [2007] ) We assumed an ini¬ 
tial pressure scale height ofh = ^/r = 0.1. This value corre¬ 
sponds to a disk temperature ofT = 100 K. We note, however, 
that the simulation of ALMA observations (see Sect. [4.1[ ) was 
based on a 3D temperature distribution that was derived from 
3D radiative transfer simulations on the basis of the density dis¬ 
tribution resulting from the SPH simulations (see Sect. [^. The 
calculations were performed using 5x10^ SPH particles. For the 
radiative transfer post-processing, we selected the distribution of 
SPH particles after 30 orbits. This is a valid compromise that re¬ 
flects the evolution of the density distribution in the inner part of 
a circumbinary disk sufficiently well. On one hand, the accretion 
rate onto the binary components saturates after 20 orbits. On the 
other, the system has lost 5% of all particles as a result of accre¬ 
tion after 40 orbits. The outer radius of the SPH simulated disks 
was set to rout = 15 AU. The total dust mass of the disk amounts 
to Mdust = 4.4 X 10“^ Mq. For investigations that have to respect 
the influence of the a complete disk, we assumed an outer disk 
following the expression in Eq. ( [Kenyon & Hartmann[[1987[ 
[Shakura & Sunyaev|1973[ ) , which was added to the numerically 
calculated inner disk. This approach is necessary to take into ac¬ 
count possible shadowing effects in close to edge-on orientations 
of the disk. 


Pout (r,z) oc 




( 1 ) 


where 


h (r) = ho 



( 2 ) 


In Eq.[^was y = 1,13 = 514, and Rq = 6a. The parameter ho was 
chosen to fit the scale height of the numerically simulated disk 
at the radius Ro (approximately ho = 0.6a). The transit between 
the numerically calculated density profile and the outer analyt¬ 
ically described region was smoothed in a range of 2a with a 
Gaussian to avoid high density gradients. Now, the outer radius 
was set to Rout = lOOAU and the total dust mass amounts to 
4^dust = 10”^ Mq or Mdust = 10”^ M©. This leads to a mass ratio 
between the disk and primary component of the binary of 0.4%, 












































Ruge et al.: Structures in circumbinary disks: Prospects for observability 



400 i 


Fig. 3. Temperature maps along vertical cuts through the density dis¬ 
tributions of the disk model with the parameters q - 0.1, i - 30°, and 
e - 0.0 shown in the top right corner in Fig.[^ 


hence effects of self-gravitation can be neglected ( |Artymowicz| 
|& Lubow|1994| ) 

Parameter space: We considered various values for the binary 
mass ratio q = ^sec/Mpnm, the eccentricity of the binary orbit, and 
its inclination to the midplane of the initial disk. For binary mass 
ratios of ^ = 1, 0.1 and 0.01, we explored binary orbits with ec¬ 
centricities of ^ = 0.0, 0.3, and 0.5. The orbits are coplanar to the 
midplane of the initial disk. Simulations of ALMA observations 
were performed when the secondary component is in apastron. 
Figure illustrates the projected density distributions of these 
disks in a face-on orientation. 

We also investigated binaries with a low-mass secondary 
component = 0.1 and 0.01) on circular orbits that are inclined 
by i = 5°, 10°, 20°, or 30° to the midplane of the initial disk. In 
these cases, we post-processed the SPH simulation output with 
the secondary component perpendicular to the line of nodes and 
underneath the circumbinary disk plane. Figure [^illustrates the 
projected density distributions of these disks in a face-on orien¬ 
tation. In addition to that, we calculated 21 snapshots of a disk 
model with ^ = 0.5, / = 0° and ^ = 0.3 to investigate structural 
changes in the disk structure during one period of the binary. 

We note that once the binary clears the inner cavity and induces 
spiral arms (e.g., q > 0.01) it is not feasible to trace from the 
density distribution back to the binary mass ratio. 


3. Thermal disk structure and the observable 
quantities of the disks 

Radiative transfer post-processing: To derive more realistic 
temperature distributions as the basis of simulated ALMA ob¬ 
servations (see Sect. |4.1| ), the heating of the disk by the binary 
radiation was simulated using the Monte Carlo 3D continuum 
radiative transfer code MC3D ( [Wolf et ar||1999| |Wolf||2003] ). 
Subsequently, reemission maps were simulated. Only, for edge- 
on orientation of the disk, we used models with an analytically 
added outer disk (see Sect. [^. 

Dust density distribution, dust properties: We used a spher¬ 
ical grid that was subdivided linearly into azimuthal (angular 
resolution: 1.9°) and polar direction (angular resolution: 3.6°). 
In the radial direction the resolution is 0.133 AU. As a result 
of the geometry of the circumbinary disks, the SPH particles 
are concentrated in 8 % of the grid close to the midplane of the 
initial disk. We substituted the distribution of SPH particles by 
a distribution of spherical dust grains that consisted of 62.5 % 
silicate and 37.5 % graphite (optical data from |Weingartner ^ 
DrainepOOT] and concept by [Mathis et al.||1977| ) with a density 



period [° ] 

Fig. 4. Fluctuations in the thermal re-emission flux of the disk during 
one period of the binary. The binary mass ratio is q = 0.5, the eccen¬ 
tricity e = 0.1. Different line types and colors represent different wave¬ 
lengths. 


of the dust material of pdust = 2.7gcm“^. Each SPH particle 
represents a dust grain size distribution follo wing the powerlaw 
d^(^dust)/d^dust ^ ^dust~^'^ ( Dohnanyi|1969t The selected range 
of grain sizes was (2dust ^ [0.005 pm, 100 pm]. Dust of this com¬ 
position and size distribution is frequently found in young stellar 
objec ts (e.g. [Wolf et al.|2003t|Madlener et al.|2012[|Grafe et al 


[201 3| ). The absolute number density of dust grains was scaled 
according to the total mass of dust in the disk. 

Central binary star: Both components of the binary system 
were considered as radiation sources. Additional possible heat¬ 
ing sources were neglected (e.g., viscous heating). The stars 
were assumed to be blackbody radiators that are characterized 
by their effective temperature and luminosity. The primary com¬ 
ponent is a Herbig Ae type star with a luminosity of L = 43 Lq, 
an effective temperature of T = 9500 K, and a stellar mass of 
M = 2.5 Mq. The binary mass ratio q defines the properties of 
the secondary component as listed in Table 


Table 1. Properties of the secondary component for the selection of 
binary mass ratios from Sect.j^ 



L 

[Lo] 

TeS 

[K] 

M 

[Mo] 

1 

43 

9500 

2.5 

0.5 

0.95 

4000 

1.25 

0.1 

0.7 

4000 

0.25 

0.01 

0.003 

2700 

0.025 


Fluctuating thermal disk structure: The asymmetries result¬ 
ing from the binary-disk interaction have affect the ability of the 
stars to illuminate the disk and therefore the temperature pro¬ 
file. The disk temperature profile is characterized by high vari¬ 
ations in vertical, azimuthal, and radial direction. We discuss 
these variations on the example of a disk model with the pa¬ 
rameters ^ = 0.1, / = 30°, and e = 0.0. Figure shows vertical 
cuts through the temperature distribution of this disk along the 
zone of the lowest and highest disk deformation ((f = 131° and 
(f = 40°, respectively). At a distance of 5 AU from the primary, 
the ratio of the surface temperatures along cuts is 2.6. The higher 
temperature along ip - 40° is the result of a better illumination 
of this area through direct stellar light. The higher temperature 
leads to a higher thermal dust re-emission in this region. 

Reemission light-curve: As a result of the variable illumina¬ 
tion described above, not only the density structure but also the 
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Fig. 5. Simulated ALMA maps at different wavelengths for eccentric 
binary orbits. The binary mass ratio is ^ = 0.1. The synthesized beam 
size of ALMA is plotted in the lower left corner of each image. The 
disks are identified by an annulus. 


temperature profile and thus the resulting reemission flux of the 
disks is time-dependent. We discuss this effect on the example 
of a disk model with a binary mass ratio of ^ = 0.5 and an ec¬ 
centricity of ^ = 0.1. For wavelengths of lOyum, 330yum, and 
1300yum, the thermal re-emission light-curve of this disk is pre¬ 
sented in Fig 1^ in 21 snapshots during one orbit of the central 
binary (1.7 yr). Although the rotation of the binary is periodic, 
there is no periodicity in the light curves at any wavelength on 
the scale of the orbital period of the binary. This behavior is con¬ 
nected to the temporal evolution of the density waves and spiral 
arms in the disk. During the orbital motion of the central stars, 
the structure of a surrounding disk changes, for example spiral 
arms are created and vanish again. Locally, this results in a vari¬ 
ation of dust mass per dust temperature interval and leads to a 
temporal dependence of the re-emission flux of the disk. The ef¬ 
fect of variations of the dust temperature on the resulting reemis¬ 
sion flux is strongest at the shortest wavelengths. In agreement 
with this, we And that the amplitude of the fluctuation decreases 
with increasing wavelength. More specflcially, the variations in 
the (sub)millimeter range are in the order of 2 %, but reach up 
to 10% at mid-infrared wavelengths. Furthermore, the magni¬ 
tude and temporal evolution of the reemission radiation depend 
on the observing wavelength, tracing regions with characteristic 
temperatures in the disk. Equivalent fluctuations in the scattered- 
light regime are observed in detail for example in the system of 
HH30 ( [Wood et al.|2000[[Tambovtseva & Grinin|2008| ) 

4. ALMA observing setup 

Since the wavelength, object position, and exposure time are al¬ 
ready flxed, the speciflc conflguration of the ALMA antennas 
is left as the only free parameter. We aim at resolving at least 
structures on an angular scale of 0.07" lOAU at a distance 
of 140 pc) and thus only consider ALMA conflgurations tat ful- 
flll this requirement (i.e., at 1.0 mm a max. baseline > 3500 m). 
Additionally, the S/N of a synthetic observation map has to be 
greater than flve for further use. Accordingly, the S/N is deter¬ 
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Fig. 6. Simulated ALMA observations for systems with binary orbits 
that are inclined with respect to the midplane of the initial disk. The 
binary mass ratio is ^ = 0.1. The synthesized beam size of ALMA is 
plotted in the lower left corner of each image. The detection of ring-like 
structures strongly depends on the disk inclination. 


mined from the background signal in the simulated observations. 
If there is more than one possibility of an observation that ful- 
fllls these criteria, we chose the simulated observation with the 
highest resolution for each wavelength. Consequently, the beam 
size is of similar size at all wavelengths. 

We simulated ALMA observations with the CASA 4.2 sim¬ 
ulator ( [Petry et ar]|20I2| ) for the (arbitrary) reference position 
of the Butterfly Star (IRAS 04302-1-2247, a = 04h33m, S = 
-r22°53', J2000) in the Taurus-Auriga star-forming complex at 
a distance of 140 pc. Thermal and phase noise awere considered. 
The detailed numerical setup for a wavelength between 740 pm 
and 1.3 mm is explained in Appendix |A| Shorter wavelengths 
were not considered because earlier studies indicated that the 
very extended ALMA conflgurations will not allow achieving a 
sufficiently high signal-to-noise ratio ( [Ruge et al.|2013| ) 


4.1. Results 

4.1.1. Circumbinary disks in a face-on orientation 

For all considered orbital elements (i < 20°) of the binary sys¬ 
tem and binary mass ratios q > 0.001, ALMA is able to detect 
the inner cavity of the disks (see Figs. and |^. In particular 
for low eccentricities and inclinations the disks appear as bright 
annuli in the simulated ALMA maps. The size of the inner cav¬ 
ity is identical for all considered wavelengths. An increase of 
the eccentricity of the binary orbits leads to asymmetries in the 
brightness distribution of the rings (Fig.[^. Spiral arms are not 
traced with this ALMA setup. The behavior of disk with bina¬ 
ries on inclined orbits (/ < 20°) is similar. Therefore it is not 
possible to distinguish between binaries on moderately eccen- 
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Fig. 7. Detectability of the spiral arm features in a circumbinary disk (model parameters: q - 0.5, i = 0°, and e = 0.3). Left: A = 740/urn, 
texp = 2 X 5h, max. baseline: ^ 15km, contour lines for signal-to-noise (S/N) level plotted in the figure. Right: corresponding midplane density 
map in arbitrary units. 


trie or inclined orbits. Only for higher inclinations (/ > 20°) is 
the structure of the disks in the simulated ALMA maps differ¬ 
ent from all other cases. As already discussed in Sect., the 
vertically warped structured of the disk leads to zones of higher 
temperatures - and therefore higher re-emission - in some of the 
outer parts of the disks. ALMA is able to trace this feature (see 
Fig. [^. For binary systems with a binary mass ratio q = 0.001 
the simulated ALMA maps do not show any evidence of a per¬ 
turbation. 


Detection of binary-induced spiral arms: In the density maps 
(Figs. andspiral arms are the most prominent feature in the 
disk. Therefore, we investigated whether it is possible to trace 
these binary-induced spiral arms with ALMA. For an example, 
we selected a disk model with the parameters ^ = 0.5, / = 0°, and 
^ = 0.3. It is characterized by a prominent spiral arm feature that 
is conserved in the re-emission maps. For the required resolution 
we focused on a wavelength of 740 jum and explored the most 
extended configurations of the ALMA array (#28 in the CASA 
array configuration files) with a longest baseline of about 15 km. 
We calculated the observing time with the ALMA Sensitivity 
Calculator, resulting in 2 x 5 h to achieve a root mean square 
noise of 0.013 mJy. 

The result is plotted and compared to the corresponding den¬ 
sity map in Fig.|7] We added contour lines to the figure that in¬ 
dicate the S/N levels. With ALMA it is possible to trace a spiral 
arm in our circumbinary disk models of 2 x 3 AU^ with the S/N 
areas that are shown in Fig.|7] 


Discussion: In general, one should be aware of tracing from 
structures in (simulated) ALMA observations back to the real 
origin of these structures. Indeed, ALMA is able to resolve and 
trace the spiral arms that are prominent in our disk models, but 
from the ALMA observation alone one will not be able to con¬ 
clude on the binary nature of the central objects. 


4.1.2. Circumbinary disks in an edge-on orientation 

Required ALMA array configuration and data analysis: Now 
we changed the disk orientation to edge-on and determined 
whether it is possible to trace features of the binary perturba¬ 
tion of the disk in the simulated ALMA maps. We focused on a 
wavelength of 1.3 mm and an ALMA configuration (#22) with a 
longest baseline of 8 km of the ALMA array. This combina¬ 
tion results in a synthesized beam of about the same size as the 
vertical extension of our disk models and maps with the highest 
possible S/N. From the simulated observations we extracted a 
horizontal cut through the disk and plot the resulting flux against 
the distance from the center (radial brightness profile). To avoid 
mistakes in the interpretation of the data that are possibly based 
on the fitting process of the analytically described outer disk onto 
the numerically simulated inner disk (see Sect. [^, we only fo¬ 
cused the radial brightness profile on a range of 20 AU around 
the center of the disk. Additionally, the radial brightness profile 
is influenced by the orientation of the disk to the observer. We 
investigated the range of inclinations of the objects from 85° to 
95° and found the following results to be valid from 88° to 92°. 
The considered dust masses of the disk are M^ust = 10”^ M© and 
^dust = 10“"^ M© (see Sect.[^. 

This mass selection marks the transiting of the disks from 
optically thin to optically thick along the line of sight from the 
observer through the edge of the disk to its center (i.e., from an 
averaged, effective optical depth: Tio-sm©, i. 3 mm = 0-65 + 0.04, 
highest value in disk midplane Tio-5Mo,max = 0-88, to an aver¬ 
aged, effective optical depth: Tio-4Mo, i.Smm = 4.16+0.50, highest 
value in the disk midplane Tio-4Mo,max = 8.79). 

The profile of the radial brightness distribution is a byproduct 
of the optical depth structure of the disk. It is thereby directly 
linked to the disk density profile. In our perfect edge-on case an 
asymmetry in the radial brightness profile therefore is the result 
of asymmetries in the density distribution of the disk, which is 
here induced by the binary. 

We find that the influences of inclined and eccentric binary 
orbits are partially different, therefore we explicitly differentiate 
between these two cases in the following. 
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Fig. 8. Edge-on brightness profile of disks around an eccentric binary. Observing wavelength: 1.3 mm. The longest baseline of the ALMA config¬ 
uration is ^ 8 km. Left: binary mass ratio ^ = 0.1; right: binary mass ratio q = 0.01. Red: M^ust = 10”^ M©, black: M^ust = 10”"^ M©. 
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Fig. 9. Midplane brightness profile of binaries seen edge-on that are inclined to the initial disk midplane. Wavelength: 1.3 mm. The longest 
baseline of the ALMA configuration is 8 km. Left: binary mass ratio q = 0.1; right: binary mass ratio q = 0.01. Red: M^ust = 10”^ M©, black: 

Mdust = 10-^M©. 


Eccentric binary orbits: For binary mass ratios q = 0.1 and 
0.01 the radial brightness profiles for eccentric binary orbits are 
plotted in Fig.[^ For ^ = 0.1 our disk models have inner cavi¬ 
ties that become noncircular with increasing eccentricity of the 
binary orbit (see Fig. [T]). The form of the cavity and the spiral 
arms, which are density waves that propagate through the disk, 
create a radial brightness profile that is not symmetric to the disk 
center. In particular, the brightness profiles of the models with 
e = 0.0 (circular cavity) and of the models with e = 0.5 (non¬ 
circular cavity) present increasing asymmetries with increasing 
eccentricity of the binary orbits. The peak flux does not de¬ 
pend on the eccentricity of the binary orbit. Additionally, disks 
with Mdust = Mq have lower peak fluxes than disks with 
^dust = 10”^ M©. This is explained by an optical depth effect 
because for the more massive disk the signal from the warm 
inner layer of the disks is hidden by the outer parts. For disk 
dust masses of Mdust = 10”^ M© the binary nature of the sys¬ 
tems leads to an increasing offset of the flux maximum with an 
increasing eccentricity. 


In the special case of systems with q = 0.01 (Fig. right) 
the binary was not able to form an inner cavity. Because of 
this, the disks with M^ust = 10”"^ M© are already optically thin 
(tio-4Mo, 1.3 mm = 0.22 + 0.04, Tio-4Me,max = 0.75). Therefore, the 
local minima in the radial brightness profile do not result from 
an optical depth effect. Their locations can be associated with 
the position of a gap formed by the secondary component in the 
density distribution of the disk (see Fig.[^. Thereby, the asym¬ 
metry of the gap is due to the eccentricity of the binary, which 
is also be traced by the radial brightness profile. Since the disk 
was optically thin before, the decrease of the re-emission flux 
for disks with dust masses of Mdust = M© can be explained 
by the decrease of the disk mass. However, the feature in the ra¬ 
dial brightness profile resulting from the gap in the disk density 
distribution is conserved. 

Since the SPH calculations are mainly sensitive to the binary 
mass ratio q, a giant planet (Mpianet = lOMjup) around a solar- 
type star (M* = 1 M©) will induce similar perturbations in the 
disk (for a detailed study see jRuge et al.|2013] ) 
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In general, the presented results are also valid for other disk 
sizes and masses as long as the optical depth along the line of 
sight from the disk center through the edge of the disk to the 
observer is similar to the values shown above. 


Inclined binary orbits: In particular, the brightness profile is 
not axis-symmetric to the rotation axis of the disk anymore (see 
Fig.|9j. In contrast to the case of eccentric binaries, the peak fiux 
increases with increasing inclination of the binary orbit for both 
considered disk masses in our setup. The reason is the warped 
surface of the disks, which leads to a b etter irradiation of some 
parts of the disk (see Sects. [3]and [m] ). 


Discussion : The radial brightness profile depends on the global 
density distributions of the disk. Thus, large-scale perturbations, 
inner holes, spiral arm features, but also the total disk mass, 
affect the shape of the profile. In addition to disk perturbation 
through binary system s, cavities can form through photoevapo¬ 
ration or planets (e .g., |Espaillat et al.||2Q12t Uribe et ar]|2011[ 
[Cossins et al.|2010| ) Therefore, it is not feasible to safely trace 
from a radial brightness profile back to its origin. 


5. Conclusion 

For a simple circumbinary disk model, we investigated the fea¬ 
sibility of tracing characteristic density structures in circumbi¬ 
nary disks that result from the binary-disk interaction. For close 
binary systems, we analyzed the impact of the binary mass ra¬ 
tio, its eccentricity, and its inclination on the appearance of the 
global disk structure. For this purpose, we simulated the 3D den¬ 
sity structure of the circumbinary disk (SPH simulations), calcu¬ 
lated the corresponding dust temperature structure and reemis¬ 
sion images (radiate transfer simulations), and derived simulated 
ALMA observations (CASA 4.2). Our main results are as fol¬ 
lows: 

- The (sub)millimeter reemission brightness distribution is 
characteristic for given binary mass ratios q. However, for 
q < 0.001 the effect of the binary nature, the resulting vari¬ 
able gravitational potential, and thus its impact, is too weak 
to create structures that are prominent enough to be traced 
with ALMA. 

- Binaries on highly inclined orbits create a warped density 
profile that leads to illumination conditions and temperature 
profiles of the disk with large variations in azimuthal direc¬ 
tion. Within the parameter space considered in this study, the 
surface temperature can locally be increased by a factor of 
2.6 as compared to unperturbed disks. 

- The time-dependent illumination leads to temporal fiux vari¬ 
ations in the thermal re-emission of the disk on the order of 
a few percent. The fluctuations are directly linked to the ap¬ 
pearance of density waves and spiral arms in the disk. 

- The inner cavity (e.g., q > 0.001) can be traced at wave¬ 
lengths > 740yum with ALMA for all disks in our setup and 
assuming standard observing conditions. 

- The simulated ALMA observations show the inner edge of 
the disks as a bright ring with local asymmetries. 

- Only at 740yum is it possible to trace the binary-created spiral 
arms of a few AU in size in the disk. 

- For disks in a perfect edge-on orientation, the simulated 
ALMA maps indicate asymmetries in the disk structure. 

We emphasize that we considered the binary-circumbinary disk 
interaction as the sole source of the studied disk structures. Other 


mechanisms, such as planet-disk interaction (|Wolf & D’Angelo 

2005 

IGonzalez et al.|2012[ |Ruge et al.||2013|), magnetic fields 

(e.g., 

Flock et al. 12014]), self-gravitation structures ([Cossins et al. 

2010 

), vortices (|Wolf & Klahr|20021), and photoevaporation (Es- 

paillat et al.||20121) may result in density waves and cavities as 


well, but were not considered in this study. The feasibility of dis¬ 
tinguishing between spiral density waves with different physical 
origins, for example through observations of different accompa¬ 
nying effects are suggested for future investigations. 
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Appendix A: CASA setup 

ALMA is a (sub)mm interferometer in the Atacama desert 
for observations in the wavelength range between 0.3 mm and 
9.3 mm. We considered the full ALMA array, which can be ar¬ 
ranged in different configurations ( [Brown et al.||2004| ) CASA 
offers a mighty toolkit for data reduction, analysis, and simu¬ 
lation for ALMA and other (sub)mm interferometers. We pre¬ 
dicted the observations with the CASA 4.2 simulator ( jPetry etj 
|al.|2012| ) To simulate an ALMA observation, CASA calculates 
the visibilities in the u-v-plane as sampled by the array con¬ 
figuration, adds noise, and finally reconstructs the image with 
the CLEAN algorithm. The observing wavelengths listed in Ta¬ 
ble [A.lj were selected to suit the atmospheric windows in the 
(sub)mm range. The bandwidth for the (continuum) ALMA ob¬ 
servations is 8 GHz, the selected exposure time is two hours. 

Table A.l. Selected wavelengths for the simulated ALMA observations. 


ALMA band 

A 

[pm] 

y 

[GHz] 

PWV 

[mm] 

Octile 

8 

740 

405 

0.658 + 0.1 

2nd 

7 

952 

315 

0.913 + 0.1 

3rd 

6 

1300 

230 

1.262 + 0.15 

4th 


Notes. PWV: precipitable water vapor 


We considered the influence of thermal noise by precipitable 
water vapor (PWV) in the atmosphere using the wavelength- 
dependent values listed in Tab. [A.l and the tsys-atm option in the 
CASA simobserve task. Additionally, we included the influence 
of atmospheric phase noise in our calculations. We assumed a 
phase-screen blowing with a constant velocity above the ALMA 


array with the PWV values and deviations listed in Table A.l 


This is numerically done in CASA with the sm.settrop task in 
screen mode. 
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